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Abstract 

We consider a class of models describing a quantum oscillator in inter- 
action with an environment. We show that models of continuous sponta- 
neous localization based on a stochastic Schrodinger equation can be de- 
rived as an approximation to purely deterministic Hamiltonian systems. 
We show an exponential decay of off-diagonal matrix elements in the en- 
ergy representation. 

1 Introduction 

The conventional quantum mechanics assumes a sharp division of the world 
into microworld of the system and macroworld of the measuring devices. Such a 
theory can be considered as a theory of an ideal measurement. However, in the 
contemporary mesoscopic physics such a division may be not adequate. We can 
imagine a direct observation of a mesoscopic body interacting with an invisible 
microparticle system. Such a record cannot be considered as a measurement 
in the conventional quantum mechanics (e.g. the behavior of the mesoscopic 
body can be reversible). Nevertheless, it gives us some information about the 
microsystem. Such a situation inspires a study of the influence of a large (possi- 
bly infinite) system on the behavior of a single particle. It seems that together 
with an investigation of a classical limit of the quantum mechanics done so 
extensively on the level of the wave function a study of the classical limit of 
the quntum theory of measurement is needed as well. It is well-known that 
we encounter difficulties with an interpretation of quantum mechanics if it is 
applied to macroscopic bodies. However, it may be that the troublesome in- 
terference principle has no observational consequences when applied to realistic 
macroscopic bodies ( apart from interesting situations on the verge of quantum 
and classical worlds Q). A research in this field ( which could in general be 
called a problem of decoherence) evoked interest recently || ||||. The prob- 
lem of decoherence seems to be closely related Bto another old problem of 
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quantum measurement theory that of the reduction of the wave packet. There 
appeared recently a conceptually simple theory of the wave packet reduction 
[|| H . In this theory the wave packet reduction is achieved through an addition 
of noise to the Schrodinger equation. Such a theory could be considered either 
as a phenomenological description of experiments or as a fundamental mod- 
ification of quantum mechanics which is supposed to describe both microworld 
and macroworld § . 

In this paper we do not intend to enter the basic epistemological problems. 
We investigate some simple models of an interaction of an oscillator with a reser- 
voir ( the simplest one is a text-book example ;see refs. H] Jl^]).The model is 
usually considered as a typical application of the quantum Langevin equation. 
We have shown in our earlier papers [flif fl5| that quantum dynamics can be 
expressed in terms of conventional diffusion processes (time-ordering replaces 
the non-commutativity). Such a stochastic description can be considered as a 
mathematical version of the Feynman path integral . We think that such an 
approach is especially useful for a description of the dynamics of a subsystem. We 
can investigate the dynamics of a subsystem through an exact elimination of the 
coordinates of the environment ( such a description is well-known for classical 
stochastic systems In principle, this is the same strategy as suggested first 

by Feynman and Vernon || and applied e.g. in refs.Q However, instead of 
using the path integral we eliminate the coordinates of the environment on the 
level of stochastic equations. We get (without any approximation) closed equa- 
tions for a particle moving in a reservoir. Our approach is based on the classical 
probability applied as a technical tool to standard quantum mechanics. It comes 
out from our study that a noise which initially enters the calculations as a tech- 
nical tool acquires properties and functions of a physical noise. The Brownian 
motion which initially is a realization of the sum over paths describing via the 
Feynman formula the unitary Schrodinger evolution acquires an interpretation 
of a dissipation when we pass to the subsystem. There is still another source 
of noise (which we call a deterministic noise). The reservoir degrees of freedom 
evolving according to the deterministic Schrodinger dynamics in the limit of an 
infinite number of degrees of freedom model quite well the white noise. This 
model of noise is similar to the mechanical model of the Brownian motion dis- 
cussed in the well-known paper of Ford et al [[t3) . We show that the Brownian 
noise together with the deterministic noise give in a proper limit the stochastic 
Schrodinger equation. As a consequence of the stochastic Schrodinger equation 
we obtain a model of a continuous spontaneous localization J^] . We show that 
in the energy representation (for a sufficiently large energy) the density matrix 
decays exponentially to its diagonal (decoherence) . A localization in momentum 
( or space) is shown to hold true only in a limit of small (resp. large ) oscillator 
frequency. 
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2 Stochastic representation of quantum dynam- 
ics 



In our earlier papers [14]- [16] we have described the unitary Hamiltonian time 
evolution of pure states in terms of the Brownian motion. This formalism can be 
considered as a mathematical version of the Feynman path-integral. On a formal 
level the method consists in the replacement of the integral over real paths q(.) 
by an integral over complex paths Viq(-) (where i is the imaginary unit). As a 
next step we change variables in the functional integral. The change of variables 
is expressed as a solution of a stochastic differential equation (depending on the 
Brownian motion). Then, the non-Gaussian functional integral can be reduced 
to the Gaussian integral over solutions of the stochastic equation. The aim of 
this formalism is to express in an explicit form the solution of the Schrodinger 
equation. Hence, it is completely equivalent to standard quantum mechanics. 

Let U t (t> 0)be a unitary Schrodinger evolution determined by the Hamil- 
tonian 

H = -h 2 — A + V (1) 
2m 

Assume we know \\ >t— Ut\x > ( x( x ) w iU denote \x > in the coordinate 
representation, where x G R n ). Let the initial condition ip for the Schrodinger 
equation be of the form ip(x) = x(x)4>(x) where tf> is an analytic function. Then, 
the solution ip t of the Schrodinger equation with the initial condition ip can be 
expressed as Xt4>t, where <f>t is the solution of the equation 

dtfa = ^A0 f + — x t _1 VXtV^ (2) 
Am m 

with the initial condition <f>. 

We express the solution of eq.(2) by a stochastic process 14 1( for time- 
dependent diffusions see ref.E^]). Then, the unitary Schrodinger evolution gen- 
erated by the Hamiltonian (1) can be expressed in the form 



(U t i/j)(x) = X t(x)E[<f>(qt(x)} (3) 

where qt{x) is a complex diffusion process starting at t=0 from x and solving 
the stochastic differential equation (here < r < t) 

ih 

dq T = —Xt-T^Xt-T{qr)dT + Xadbr (4) 



where 



and 



h 

a = 

m 
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The Brownian motion bt is defined as the Gaussian process with independent 
increments and the variance 

E[b 2 t ] = t 

In order to express the solution of the Schrodinger equation for negative time 
we may apply the complex conjugation 

Ut$ = 1>-t 

The action of arbitrary operators on states as well as correlation functions 
in arbitrary states can now be expressed as expectation values with respect to 
the Brownian motion |l4|| .For example 

< X \F(x t )G(x)\ X >=< x\UtF(x)U t G{x)\ X >= E[ f dx\ X t(x)\ 2 F(x)G(q t (x))} 



(5) 

When F = 1 then eq.(5) expresses a generalization of the notion of an invariant 
measure to time-inhomogenous processes ( if x is a stationary state then |%(:r)| 2 
is exactly the invariant measure for q t fill ) i.e. 

E[J dx\xt(x)\ 2 G(q t (x))} = J dx\ X (x)\ 2 G(x) 

Let us consider as an example a model of independent oscillators with m = 1 
(here and later we do not specify the number of oscillators, it can be infinite) 

fc K k 

The ground state solution of the Schrodinger equation with the Hamiltonian (6) 
reads 

Xo{x) = exp(-—J2uj k x 2 k ) (7) 

k 

The stochastic equation (4) takes the form 

dqtc = -iuj k q k dt + Xadb k (8) 
The solution of eq.(8) with the initial condition x reads 

q k (s) = exp{-iuj k s)x k + Act / exp(-iuj k (s - r))d& fe (r) (9) 
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3 An interaction with the environment 



The formalism of sec. 2 can easily be applied to Hamiltonians with explicitly 
known x*-An example of such a class of models are Hamiltonians ( often en- 
countered in quantum optics) of the form 

H S = Y J KM,A + )A+A U (10) 

where is the anihilation operator,yl+ the creation operator and the Greek 
index indicates the space coordinate. In the case of the Hamiltonian (10) the 
harmonic oscillator's ground state (7) is also the ground state of Hg. Hence, we 
may take \o as Xt in eq.(3). 

We add to Hg the coupling to the reservoir. In this paper we restrict our- 
selves to the simplest one-dimensional model corresponding to h = 1. In such 
a case we obtain a standard model of an interaction with the environment (see 
[|I),|]l2) and Q for a recent review). 

H = HluoA+A + hu k a+a k + ^g k A + a k + gla+ A (11) 

k fc 

Using the Schrodinger representation of the creation and anihilation operators 




we express H in the form (1). It is clear that the product of oscillator ground 
states is the ground state for H.In order to simplify the discussion we assume 
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that 

9k = ivk 

are purely imaginary ( if g k had a real part this would change only the form of 
the noise in the stochastic equation (4)). Then, taking as \t in e Q-(4) the ground 
state (of Q and q k oscillators) we obtain the following equations 

dq k = -iuo k q k dt + v k . —Qdt + Xadb k (12) 
V Wfc 

El ^k 
v k * — q k dt + Xadb (13) 

k v w ° 

We solve first eq.(12) for q k . We obtain 



q k (t) = exp(-iu) k t)x k + v k ^J^- J Q exp{-iu k {t - s))Q{s)ds+ 
+Acr J Q exp(-iuj k (t - s))db k (s) 

Let us introduce the notation 



(14) 



n d{s) = v k x k . —exp(-iu k s) (15) 



and 



/ — i — r s 

N R (s) = Act VVw— / exp(-iu k (s - T))db k {r) 

u V ^0 Jo 



(16) 



corresponding to a "deterministic" and random noise. We are interested in a 
computation of correlation functions in eq.(5) with \t defined in eq.(7). Then, in 
such computations x k play the role of independent Gaussian random variables 
with the covariance 

E[x k x r ] = - — S kr (17) 
2uj k 

Hence, the noise No is a Gaussian complex stochastic process with the covariance 

E[N d (s)N d (t)} =-£-J2vlexp{-iu k (T + s)) (18) 
2uJ ° k 

E[N d (s)N d (t)} = JLy^tftexpi-iukiT - a)) 

This noise is of the same order as the Nr noise (16). Next, let us denote 

K(s) = ^v 2 k exp(-iu} k s) (19) 
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We insert the solution qt of eq.(14) into eq.(13) for Q. We obtain an equation 
describing an oscillator moving in a reservoir 

dQ = -iuoQdr - [ K(t - s)Q(s)dsdT - N D {r)dT - N R (r)dT + Xadb (20) 
Jo 

It can be seen from eq.(19) that K(s) is positive for small s, hence it fulfills the 
role of a dissipation. Np ,Nr and b add to a random external force. 

4 Stochastic Schrodinger equation from an in- 
teraction with the environment 

In refs.[l]-[3] a classical behaviour of quantum systems of the type (11) (in par- 
ticular the decoherence) is explained by means of the Feynman- Vernon influence 
functional. This functional results from an average over environmental degrees 
of freedom. We could derive it explicitely from the solution of eq.(20). There is 
another theory of decoherence. The authors |^[^) explain the emergence of clas- 
sical properties through an addition of some stochastic terms to the Schrodinger 
equation. In spite of the ardent opposition of the two groups ( see the letters 
in April 1993 issue of Physics Today |2lJ]) we would like to reconcile the two 
approaches. We restrict ourselves here to the model (11). In eq.(20) the time 
evolution of the coordinate Q is influenced by the environment through the 
additional noise terms Nr and Njj. Nr (eq.(16)) can be considered as a ran- 
domness resulting from the Feynman's sum over paths jl6| ( weighted sum over 
paths instead of the classical extremal path). The "deterministic" noise No ( 
eq.(15) ) is produced by the environment of oscillators evolving in a determinis- 
tic way according to the Schrodinger evolution determined by the Hamiltonian 
(11). We can show that similarly as in the model of Ford et al (llj the environ- 
ment can influence the system as a random force of the form of the white noise. 
We note first that if in eq.(18) Vk «const and ~ k then the noise Nd(s) is 
strongly correlated only at coinciding arguments. Under these assumptions on 
Vk and iOk the dissipation kernel K(s) in eq.(19) can now be approximated by 
(such a behaviour of K is sometimes called the Ohm dissipation) 

K(s) = aS(s) 

where a is a positive constant. In fact, we get such a formula if u>k — few, where 
k € [0,oo], and vi — -ujAk. In such a case the sum over infinitesimal Ak can 
be replaced by an integral 

dkexp(-ikujs) = —S(s) 

Then, the term s _1 is considered as small in comparison to S(s). The cor- 
responding choice of v\ leads to the following approximation to Njj (again the 
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term (t — s) 1 is neglected) 

E[N D (s)N D (t)} = (21) 



E[N D (s)N D (t)} =heS{t-s) 

where 

a 
2uj 

So,Nd is the white noise. We denote the corresponding Brownian motion by 
B D i.e. 

as 

Next, let us consider the "quantum" noise N R (eq.(16)) in more detail. Let us 
compute its correlation functions 

E[N r (s)N r (t)} = A V v 2 k {exp{-iLo k \T - s\) - ex V (-iuj k (j + a))) (22) 
2lj - 

£[-/Vr(s) A^(t)] = hmin(s, t) — V] u k v 2 k exp(iuj k (s - t)) 

^ I 

In the approximation of the Ohm dissipation we get by the same arguments as 
applied at the derivation of eq.(21) 

E[N r (s)N r (t)] - hed(s - r) (23) 

Eq.(23) means that N R is the real white noise (the second of the expectation 
values in eq.(22) does not enter any subsequent formulas). We denote the inte- 
gral of N R by B R (the real Brownian motion independent of Bp). Now,eq.(20) 
reads 

dQ = -ituoQds - aQds - dB D (s) - dB R {s) + Xadb(s) (24) 

( all the noise terms on the right hand site of eq.(24) are independent). 

The dynamics of a subsystem is usually investigated in terms of the time 
evolution of its density matrix. We define the density matrix of a subsystem ( 
with coordinates X ) as 

p P {X,X') = J]Jdx kP (X,x;x,X') (25) 

k 

We make the usual assumption that initially the system and the environment 
are separated i.e. p is of the product form 



p(X, x- x, X') - exp(-^(X 2 + X' 2 )MXMX')po(x, x) 
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We have expressed the initial wave function of the subsystem as a product of 
the ground state wave function and arbitrary v{X) ( this is a proper choice 
for the class of models (10)-(11)). We choose for po the probability density 
in the ground state (7). This is a reasonable assumption for the reservoir at 
zero temperature ( the stochastic formalism of sec. 2 at positive temperature is 
discussed in 20 ). 

We can express in a simple way the evolution of subsystem's density matrix 
in terms of the Brownian motion 

p P (t;X,X')= 

/n fe ^ea ;p (-^(x 2 + x' 2 ))| Xo (x)| 2 i ! ;KQ t W)]i?[KQ t (x'))] 1 ' 

For the model (11) the stochastic process Q is the solution of the equation (20). 
When we apply the definition of the noise Np> (eqs.(15) and (17)) then the 
integration over Xk can be expressed as an average over Np> 



pp(t;X,X')^ e ^(-^(X 2 + X' 2 ))ii; D [i?[KQt(X))] J B[KQ i (X'))]] (27) 

where E[..] means the expectation value with respect to the Brownian motion b 
and Ep, means that the integration over x is replaced by an expectation value 
with respect to the complex white noise Np>. An application of the standard 
Ito stochastic calculus ( see e.g.|l7|]) shows that v{Qt) fulfills the Ito equation 



dv = -^(—iujoQds — aQds — dBp>(s) — dBp + Xadb)- 



(28) 



The wave function ip(X) = exp(— ^^—)v(X) satisfies the stochastic Schrodinger 
equation 

# = -~H Q il)ds - ^H ipds - Li)(dB D (s) + dB R (s)) + XoL^dbls) (29) 
n n 

where Hq is the Hamiltonian of the oscillator (the formula (6) for the luq- oscil- 
lator) 

T - 9 

T)X ~h 

We can derive the master equation for pp from eq.(27) without any further 
approximation. For this purpose we differentiate pp in eq.(27) and apply eq.(29). 
Then, elementary rules of the stochastic calculus JllJ lead to the formula 

dtpp = -~[i?o,pp] ~ ^{H pp + ppH ) + heL PP L+ (30) 
We have got the master equation exactly in the Lindblad form |2^| as 

h 2 , 

Ho = —L+L 



The approximation which leads to the stochastic equation (29) is a Markovian 
approximation implying linear differential equations for the density matrix.The 
Lindblad form of the evolution equation ensures that the evolution preserves 
the normalization of the density matrix and its positivity. 

We consider eq.(30) in the energy representation. Let \k > be the k-th 
eigenstate of H . Then, eq.(30) for matrix elements reads (for a derivation it is 
sufficient to notice that L is proportional to the anihilation operator) 

d t < j\p p\k >= -iiVojj - k) < j\p P \k > -eu (j + k) < j\p P \k > + 
+2e^o V / (.7' + l)(fc + 1) < 3 + MPp\k + 1 > 1 ' 

For large j and k and a wide class of operators pp (in particular for operators 
whose matrix elements have a limit ft, — > 0) < j\pp\k >~< j + l\pp\k + 1 >. 
Applying this approximation to the last term in cq.(31) we obtain the solution 

< j\pp(t)\k >~ exp{- l -{Ej -E k )t- \{\f^i- VE~k?t) < j\p P (0)\k > (32) 

here we denoted by E k the k-th energy eigenvalue. Eq.(32) shows that for macro- 
scopic energies the density matrix tends exponentially fast (with the increase of 
time) to its diagonal. So, if we take in eq.(26) v = X)p*lp c p^p> (where \p >= 7i p 
is the Hermite polynomial of order p) then for large P and large t the density 
matrix tends to Y^Lp \ c p\ 2 \p >< p\- We can see that the interference disap- 
pears on a time scale proportional to the inverse of the energy difference of the 
eigenstates in the sum v. Hence, there is no interference of states with macro- 
scopically distinguishable energies. There is still another interpretation of the 
result. Assume that in the framework of quantum mechanics of closed systems 
with the conventional wave packet reduction postulate an apparatus performed 
a measurement of the energy eigenvalues, but the results are unknown. Then 
,the state of the system after the measurement is described by the solution (32) 
at large time. In this sense the interaction with an environment produces the 
same effect as a measurement. 

If we consider eq.(30) in the momentum representation and assume oj$X 2 ~ 
,so that we can neglect the oscillator potential energy in comparison with the 
kinetic energy, then from eq.(30) we get approximately 

d t < p\p P \p' >~ ~{p - p'f < p\p P \p' > (33) 

Eq.(33) means a localization in momentum. If on the other hand we assume 
that ujqX 2 is large in comparison with the kinetic energy, then neglecting the 
momenta in eq.(30) we get in the position representation 

d t < X\ PP \X' >~ ~^(A - X') 2 < X\ PP \X' > (34) 
i.e. a decoherence in space. 
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5 Summary and discussion 



We have discussed a description of the dynamics of a subsystem which clearly 
displays the dissipation and noise in the subsystem. It is shown in a simple 
model that some classical properties can emerge in quantum systems from an 
interaction with a large (quantum) environment. In particular, we have derived 
a stochastic Schrodinger equation discussed in the theory of measurement. 

The theory of an interaction of microscopic and macroscpic bodies based 
on the stochastic Schrodinger equation has some appealing features. It de- 
scribes the dynamics of decoherence and localization in a model-independent 
way.However, we think that the microworld and macroworld will meet in the 
conventional quantum mechanics. In such a case the problem of the consis- 
tency of the deterministic ( in the sense of the Schrodinger time evolution) and 
stochastic evolutions will emerge. This was our motivation for an investiga- 
tion of the stochastic Schrodinger equation as a consequence of a deterministic 
model. The stochastic Schrodinger equation derived in sec. 4 seems to be a 
realistic approximation to the time evolution in a reservoir. Surprisingly, it 
resembles a model of ref. j^] for a continuous spontaneous localization during 
the energy measurement. A measurement of another observable needs another 
environment which in general would be difficult to define explicitly. However, it 
is important to enquire whether we can reduce the number of basic principles of 
quantum mechanics without any modification of the Schrodinger equation and 
the superposition principle. 
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